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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 

ELECTRIC  CONDUCTIVITY  OF  WEAKLY  IONIZED  GASES 

In  calculating  the  electric  conductivity  of  gases  and 
plasmas,  it  is  not  sufficient  to  define  only  the  effective 
collision  frequency.  Therefore,  two  independent  quantities 
are  defined  which  permit  simple  adaptation  of  existing 
solutions  of  Maxwell  equations  (at  constant  collision  fre- 
quency) to  any  realistic  plasma  in  the  presence  of  static 
magnetic  fields.  A detailed  review  of  the  theory  of  elec- 
tric conductivity  of  weakly  ionized  gases  herein  is  presented 
and  supplemented  with  many  additional  derivations  and  compu- 
tations. Included  are  tables  of  the  effective  quantities 
needed  for  the  electric  conductivity  for  an  idealized  single- 
component weakly  ionized  gas  whose  collision  frequency  has 
a power  law  dependence  on  the  velocity  (with  the  exponent  in 
the  range  of  -3  to  +3).  The  average  electron-neutral  collision 
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INTRODUCTION 


aquation  can  be  found 


-3£r 


ui/v,  ibp/v  or  Up/®,  Co  be  valid  if  Chose  si 


Florida  on  Microwave  Propagation  Through  Plasma.  Various  tables 


l.  MEAN  VELOCITY  AND  COLLISION  FREQUENCY  OF 

ELECTRONS  ( LORENTZ  THEORY) 

In  lectures  at  Columbia  University,  the  great  Dutch 
physicist  H.A.  Lorentz1  considered  the  motion  of  bound  electrons 
under  the  influence  of  an  applied  alternating  electric  field  and 
obtained  an  expression  for  the  mean  value  of  displacement.  We  will 
derive,  in  an  analogous  way,  the  mean  value  of  velocity  for  free 
electrons  under  the  action  of  electromagnetic  waves  passing  through 
plasma.  In  the  absence  of  an  applied  external  field,  between  colli- 
sions, the  motion  of  j-th  electron  n the  electromagnetic  field  is  de- 
scribed by  the  differential  equation  (Lorentz'  force  equation) 

m.  |5i  = eg,*-1"*  t f (V  * B)  i v,  = ||i  , (2) 

where  is  the  mass  of  an  electron,  E is  the  electric  vector, 
which  is  assumed  to  be  varying  harmonically  with  time  (e.g.,  micro- 
waves), £*  = Gce  *‘bC,  and  where  the  term  ^(v.  * HJrepresents  the 
total  magnetic  force  on  the  j-th  electron.  The  contribution  of 
moving  ions  and  electrons  into  the  ff  we  assume  as  neglibible. 

If  is  an  external  static  magnetic  field  then  by  solving 
Eq.  (1),  we  will  see  from  Eqs.  (94)  - (96)  thac  we  can  drop  the 


magnetic  force 


I - e - , as  long  as  Che  electron-cycloCron  frequency  i 
the  frequency  of  the  electromagnetic  wave. 

The  solution  or  Eq.  (1),  neglecting  the  las 
Vj  = ||l  = e"1"*  + C,  = ^ (e'ialt  - e'1' 

The  integration  constant  Ci  is  so  determine 
t] , immediately  after  its  latest  collision,  the  meai 
aged  over  all  electrons  j = 1 to  N, ) is  zero.  To  ti 
random  (thermal)  velocity  is  superimposed  the  order 
component  of  the  velocity  in  the  direction  of  t 
field  E.  Other  particles  (neutrals  and  to  a good  approximation  t 
heavy  ions)  are  moving  at  thermal  velocities.  Since  we  assume 
vthermal  v an<*  at  t = t,  all  directions  of  the  displace- 
ment and  velocity  are  equally  probable  (isotropic  scattering), 
then  in  order  to  obtain  mean  values  of  v we  require  v(t,  ) : 0. 


> the  electron's 
lered  (oscillatory) 
! applied  electric 


= l£!s.  e"1®1  £ [l-e*1®0*] 


This  is  the  mean  value  of  the  oscillatory  velocity,  at  a 
definite  instant,  t.  How  at  time  9 let  N, ( 6)  of  the  electrons 
still  be  moving  without  having  had  a collision.  If  v is  Che  velo- 
city independent  collision  frequency  of  electrons  wich  other 
particles  and  ' - 1/v  is  the  average  time  between  collisions,  then 


ie/T  is  the  fraction  (probability)  that  collision  will  occur  in 
the  time  interval  A0(A0  < t ) , and  N,A0/t  is  the  number  of  electrons 
that  will  (has  probability  of)  collide  and  drop  out  of  the  remaining 
group  of  electrons.  Therefore,  S',  is  changed  by 

dN.  = - N.vde  , (6) 

which  integrates  to  the  number  of  electrons  remaining  in  the  group 


The  number  of  electrons  that  will  collide  between  time  0 and  0 de 
is  thus  given  by 


The  fraction  from  the  freely  traveling  electrons  (from  Che 
original  selected  N,  ) which  will  collide  in  the  time  between  0 and 


(9) 


which  is  also  the  probability  that  the  time  before  the  next  colli- 
sion will  occur  in  the  lime  do,  and  also  the  probability  that  any 
electron  still  has  time  0 to  travel  before  the  next  collision  in 
the  time  0 and  0 + d0. 


The  average  velocity  among  the  freely  traveling  fraction 
of  originally  selected  M,  electrons  in  the  time  between  0 and 


dv  = {V>N,V  e'V0d9  • (10) 
hs  the  time  progresses  and  more  electrons  drop  out  of  the  selected 
group,  the  average  values  will  approach  (summed  over  all  N, 


electrons  freely  traveling  from 


v=  f (?)Nive‘v0dB  . (11) 

We  could  do  the  calculations  at  a different  time  ta  and  with  a 
different  number  of  Ng  electrons  colliding  at  this  time,  and 
follow  in  time  those  out  of  this  group  which  do  not  collide,  but 
since  our  solutions  to  Eqs.  (6)  and  (8)  are  independent  of  the  num- 

deed  the  proper  average  quantities.  The  mean  value  of  the  oscilla- 
tory (drift)  velocity  of  an  electron  is 


This  mean  velocity  of  a free  electron  was  derived  previously  by 


Appleton  and  Chapman^  who  followed  closely  the  derivation  by  horentz, 
and  has  been  discussed  more  recently  by  others.  ' For  a velocity- 
independent  collision  frequency  and  a weak  field  theory  (Eg  small), 
one  can  obtain  the  same  result  from  the  Boltzmann  equation,  i.e., 
neglecting  terms  of  order  Eg3 , etc?”13 

One  can  now  ask  which  differential  equation  will  give  the 
same  average  result,  i.e.,  what  form  has  the  differential  equation, 


m.  ^4m.v(?j  = eE0e-i»t  . (15) 


2.  TIME  AVERAGES 


Under  < > we  will  understand  average  quantities  over  time  t, 
t>^.  We  have  (sin2  wt>  = (cos2  <ut)  = T J (sin2ujt)  dt  ~ 1/2, 

(cos  tot)  = (sin  Hit)  = 0.  Thus  the  root  mean  square  velocity  is 
from  Eq.  (15), 

u = V < (ft?  > = ^ ° - = ^=-  Hmmr-  <16> 

The  average  kinetic  energy  imparted  to  the  electron  from 
the  field  (drift  energy)  is 

U = (Re  (eE)  • Re  (v))/g  = <^Eq  cos  wt)|v|>  , (17) 

which,  after  substitution  of  v from  Eq.  (0),  gives 

U = I m °-a)  = me  u“  = ^ me  |v|ama)[  . (is ) 

This  can  be  shown  also  from  the  Boltzmann  equation  (see  Appendix 
in  Reference  60). 


3.  PLASMA  FREQUENCY 


For  an  infinite  homogeneous  plasma  at  rest,  let  us  momen- 
tarily displace  N electrons  per  cma  along  some  configuration.^ 

A restoring  electric  field  E will  be  set  up,  which  can  be  obtained 
from  one  of  Maxwell's  equations 


Just  displacing  electrons  the  same  distance  along  one  dimension  in 
an  infinite  homogeneous  slab  (the  ions  are  assumed  to  be  stationary) 


(20) 


where  Ex  is  the  electric  field  arising  from  the  displacement  of 
electrons.  Displacing  the  same  distance  radially  N electrons  in 
an  infinite,  homogeneous  plasma  along  an  infinite  cylinder  in  a 
radial  direction,  and  using  cylindrical  coordinates 


(21) 


Displacing  N electrons 
the  radius  of  a sphere  a 


n an  infinite  homogeneous  plasma  along 
d using  spherical  coordinates  we  get 


a solution 


ie  equation  of  motion  of  an  electron  (neglecting  collisions) 
ter  the  displacement  to  the  distance  ?m-x  at  time  t = 0 
15 ,16  ^ pj3Sma  electron  oscillation) 


dt=  ° 

a solution 


The  (angular)  plasma  frequency^  o 


electron  is  given 


12 


* 10“ 


= 2nf,  , 

f = 8.978  x 10s  \rN  cycles/sec  , (27) 


where  N stands  for  the  number  of  electrons  in  1 cm3 , m is  the 
mass  of  the  electron,  e is  the  charge  on  the  electron.  Thus  the 
electrons,  neglecting  collisions,  will  continue  to  oscillate  with 


an  angular  frequency  ut/./a~  along  an  equilibrium  point;  the  slab 

sphere  with  is, / ,/T.  ^ Taking  collisions  into  account,  the  result- 
ing oscillations  will  be  damped,  transferring  some  of  the  electro- 
magnetic energy  into  energy  of  plasma  (temperature  increase, 
atomic  excitation,  etc.). 

Often  in  the  literature  plasma  frequency  is  referred  to 
as  Langmuir  frequency,  since  Tonks  and  Langmuir  were  first  to 


introduce 


4.  THE  ORDINARY  ELECTRIC  CONDUCTIVITY  OF  PLASMA 
IN  THE  ABSENCE  OF  MAGNETIC  FIELDS 

The  Lorentz17  formula  for  the  (complex)  electric  conduc- 
tivity in  the  absence  of  an  external  magnetic  field  is  deduced 
from  Eq.  (12)  and  from  the  definition  of  o,  i.e.,  from 

where  j is  the  current  density  and  R 6(3")  is  the  real  part  of  j. 
J = oE  = N,e(v) 


ar.d  wtere  w is  the  incident  angular  frequency  of  electromagnetic 
wave,  w = 2nf,  v is  the  so-called  average  collision  frequency, 


i.e.,  average  number  of  collisions  which  an  electron  makes  with 
other  particles,  or  and  cj  are  the  real  and  imaginary  parts  of 
complex  conductivity,  respectively.  (We  are  neglecting  here 
collision  of  ions  with  ions,  electrons  with  electrons,  motion 
of  ions  and  dependence  of  v on  electron  velocity.) 

The  specific  impedance  is  usually  defined  as 


is  called  the  resistivity  and 


Using  the  electric  charge  in  electrostatic  units,  me  In 

Conductivity  in  other  units  are  related  as  follows, 

S o(  sec'1 ) = 9xl09  o(mho/m)  = 9x10 20  o(abmho/cm) 


o(mho/cm)  - 


a( ahmho/cm) . 


n ? and  R where  J is  che 


«/+»e^=eto0 os»t 

= ^C,e°B“t  + "Sln“ 


wlth  <^^.TiS»tiea 

mt+  me^= 
♦-*£ 

” me(»W) 

? = Re(aE) 


Vc^ 


V*2-2) 


teu't' = "e  ■ 

ESJot- Thecorrei 

?•£  = (Re  3)  . (Re  S) 


16 


conductivity,  T|  = 1/0  for  the  direct  current  (DC)  conductivity 
(conductivity  for  u)  = 0).  The  Lorentz  formula  gives,  in  this  case 


defined  by  the  relation 


where  (v)  is  the  electron  drift  velocity.  From  Eq.  (8), 


(33) 


5.  THE  ELECTRIC  CONDUCTIVITY  OF  WEAKLY  IONIZED  GASES 
IN  THE  ABSENCE  OF  EXTERNAL  MAGNETIC  FIELDS 
(WEAK  ELECTRIC  FIELDS) 

(MAJUMDAR  - MARGENAU  THEORY) 


Ma  j uirdar  8 


the 


itivity  (or  actually  mobility)  of  plasma 


later  independently  by  Margenau19  for  the  electric  conductivity 


dv  = 4nv*dv 


v(v)  , 


j with  o give 


> etc- and  "A»'  «»* a 


■=vfeff  :-r-‘- 

EvT££!?“  ”“"4”) «“  si*)  ‘4d“, 


Introducing 


The  real  and  imaginary  parts  of  the  dielectric  constant  a 
therefore  given  by 

er  = 1 - TT  = 1 - -T=r  f ~ , W-Jw  e'W  dW 

J'A'  p J up  + Vs 

e = o^u  - -t-  „ • f -^L  e'W  W-/vT  dW  , 

1 r iJF1  p J „P  + Vs 


from  Eg.  (35)  and 


If  v is  independent  of  velocity 
| e-W  dW  = ijOT 


ordinary  Lorentz-Appleton-Champman-Langevin  relation. 


Furthermore,  the  complex  dielectric  constant  is 


(42) 


ie  real  dieleatric  c 


and  the  imaginary  dielectric  constant  is 


For  the  Maxwellian  distribution  In  the  theory  for  weakly 
plasma,  from  Eq.  (36)  and  Eq.  (33),  v=  v(W), 


= -2— 7 [ 

V*  J 


ionized 


and 


ITdLZiZ  °Z 

n frequency*  and 

- meva/2kT J dv 

" W JQ  V<W)  ^ 6 ^ ’ 

-»■>  . 

In  Eqs.  (52)  and  (53)  the  defined  effect 

ye  quantities  bee 

Mu-f 


i.e.,  a,  - 1,  a,  - 1,  a.  - 1,  as  w/v  - 

lissi~= 


23 


(50) 


and  if  u)  = 0 (direct  current  case)  from  Eq.  (33)  and  Eq.  (36), 


-i-  = -L-  f " i W4'  3 e'W  dl 

V s^rj  v 


-ic  " *F  (”^)  1 


We  will  see  later  that  for  a highly  ionized  plasma  Eos.  (50) 
and  (51)  do  not  hold  due  to  the  neglect  of  electron-electron 
collisions. 

For  uj  » v we  can  define  an  average  effective  collision 
frequency  v,  and  for  ui  = 0 an  average  effective  collision  fre- 
quency, \^jC.  For  any  other  range  of  frequencies  (e.g.. 

10““  v < ui  < 10a  v),  there  is  no  meaning  to  the  average  collis- 
ion frequency.  (It  is  useful  only  as  a theoretically  computed 
convenient  normalization  parameter  in  the  equation.)  Xn  that 
case  there  are  only  two  meaningful  parameters:  e and  e,  (or 


o and  c • , e and  a.,  or  e.  and  o ),  or  for  convenience  only 
di  and  os  (or  Hi  and  Oo , or  oa  and  a^)  defined  in  the  next 
paragraph.  From  an  experiment  in  microwave  propagation  in  a 
homogeneous  plasma  one  can  deduce  e's  (or  o's)  and  only  for 
m » v and  ui  « v compare  directly  with  the  theoretical  ~v 
or  . For  intermediate  frequencies,  using  theoretical  e's 
(or  o's,  or  a's),  one  can  also  deduce  from  experimental  e's 
(or  o's)  the  semi-empirical  average  collision  frequency  ~ 
for  the  case  of  o>  » v and  compare  this  with  the  experimentally 
obtained  or  theoretically  computed  value  at  u>  » v. 

From  Eq.  (S4)and  by  analogy  to  Eq.(48)for  any  distribution 
function,  for  u » v,  v is  defined  by 


(51a) 


§■  f 


(4,=?» 

{*)-■>•$  ■ 


Let  us  define  the  following  quantities,  which  we  will  use  frequently 


The  dielectric  constant  becomes  (see  Eqs.  ( 12),  (52)  and  (53)) 


and  the  electric  conductivity 


The  v is  so  chosen  that  o.,  a-  and  approach  unity  for  an  infinite 
w/v.  In  this  way,  o\,  a.,  and  a,  indicate  departure  from  Lorents' 
ordinary  conductivity.  For  Maxwellian  distribution ,we  already  defined 


27 


1 


(61) 


which  will  be  shown  later  is  a good  expression  not  only  for  a weakly 
ionised  gas  but  also  for  a partially  and  fully  ionized  plasma  (as 
long  as  w » v)  and  a convenient  normalisation  parameter  for  other 


tity  v, / f i but  obviously  it  cannot  be  uniquely  defined  to  insure 
compatibility,  except  only  for  m » v and  for  la  « v- 


in  the  electric  conductivity  tensor,  when  the  plasma  is 
in  a static  magnetic  field(Eq.  (36),  etc.). 

The  value  of  the  here  introduced,  and  later  tabulated  (Table 
»-l.  Figs.  C-l-C-9 ) effective  corrections  (t^.Oj.Oj)  is  that  these  allow 
us  using  any  previously  solved  cases  of  propagation  of  electromag- 
netic waves  through  a given  plasma  by  directly  interpreting  the 
used  m,/u!  and  u/v  as  actually  effective  quantities.  Thus,  if  some- 
body previously  solved  numerically  homogeneous  plasma  slab  using 
Lorentz'  conductivity  for  mp/m  = 1.0  and  gj/v  = 0.09  (with  constant 
v)  S we  can  still  use  the  results  for  velocity  dependent  v by  in- 
terpreting the  parameters  as  (mt/m). ,,  = 1.0  and  (ti>/v),ff  = 0.09. 

Then,  for  example,  for  a particular  plasma  (e.g.,  v ~ v^,  Table 
A-l)  for  a given  m,/ui  = 1/0.6741  one  obtains  o,  = 0.6741  and 
ui/v  = 0.025.  Thus,  the  reflection,  transmission,  etc.  of  the  plasma 
having  v ~ v and  ai,/ai  = 1/0-6741  and  ui/v  = 0.025  are  exactly  the 
same  as  for  a plasma  having  v = constant,  oit/iu  = 1 and  m/v  = 0.09. 


frequencies  (Ginzburg6  calls  this  v,  v,f(;  his  Eq.  (69)). 


experimental 


previous  section  from 
for  all  u!  (except  for 


iu  » v and  « = 0)  one  can  with  the  help  of  theoretical  a's  deduce 
v for  0)  » V.  There  is  no  significance  to  the  average  collision  fre- 
quency at  intermediate  s . One  can  compare  the  measured  ratio 

(as  interpreted  from  microwave  measurement  for  an  assumed  geometry 


with  the  theoretical  one. 


where  1 ^ is  the  collision  frequency  obtained  assuming  Lorentz  theory, 
or  from  the  theoretical  or,  one  can  deduce  semi-empirical  v at 
u]  » y.  Note  that  the  above  ratio  is  independent  of  electron  con- 
centration. The  other  comparison,  of  course,  is  that  of  €,  and 
6,  directly  with  the  computed  €,  and  6,  (e.g.,  through  the  use  of 
tables  of  o's).  From  experimental  €r  and  €,  in  a homogeneous  plasma 
one  can  deduce  electron  concentration  from  theoretical  of, , i.e., 
from  Eq.  (59)  taking  the  ratio  of  1 - to  and  substituting  for 

= -jpr  in’  gr ^ . (63) 


electron  concentration  obtained  assuming  Lorentz1 


= 1).  Eqa.(62)  and  <63)v 


by  F.  Te 


comparing  his  experimental  results  obtained  in  a shock  tube  for 
high  temperature  air  with  our  calculations  (see  Figs.  C-7  to  C-9). 


- ^+'’,5  V 1 'V  J}  K,  f64) 

one  obtains 

0' = W^=^r?K«  ’ 

where  v could  be  any  constant  collision  frequency.  Explicitly , for 
a weakly  ionized  homogeneous  plasma  and  Maxwellian  distribution 

where  = Ka(u>/v)  for  a given  v = v(v)  and  similarly 

K'  = Iff  / "-/|±4-  e-V/2dW  = Ke(m/v)  . (66) 

For  v = v,  = K = 1.  Furthermore,  the  previously  defined  I1  and 


th  are  independent  of  the  choice  of  v.  The  same  holds  for  ou 
Ginzburg0  quotes  and  X only  for  a few  selected  cases. 


tabulated  (for  selected  cases)  g 


Shkarof  sky2'  defined  and 


°'  ■ 

The  various  correction  factors  are  related  to  each 
the  following  relations* 


"l  ' g;  " IT" 


(69  ) 


(70) 

(71) 

(72) 


(73) 


the 


(76) 


for  fully  ionized  plasma  is  given  by  Spitzer?6  The  v 


presence  of  stationary  magnetic 
electron  cyclotron  frequency. 


by  us  extensively  for  weakly 


ionized  plasmas 

in  appendices  D-l  and  E and  for  fully  ionized  plasmas  in  Refs. 


60,  61  and  62. 

Molmud35  defined  a complex  quantity  g = gr  + ig,  (in  sec-1 ) 
so  that  one  could  formally  use  Lorentz  (Langevin)  equation  (see  Eq. 
(«), 


irHJr 


(77) 


g = + i»[i  - d,"3]  (78) 

For  in  »v,al=i%  = 1/3=1  implies  again  the  ordinary  Lorentz 
equation  with  g = v.  For  to  = 0 we  obtain  Lorentz  equation  with 

Whitmer  and  Hermann  have  lately63  introduced  vfiff  and  (o>p)eff> 
(for  weakly  ionized  gases  only),  which  are  related  to  our  published 
before31’32  cr's  (for  any  plasma) 


(80) 


These  cannot  be  easily  extended  in  the  presence  of  magnetic  fields 
as  our  a's  can  (see  Eqs. (143)-(146))  and  add  by  the  name,  we  feel, more 

frequencies  (see  Eq.(49a)  and  Eq.  (185)). 


ARGAND  TYPE  DIAGRAM 


COMPLEX  PIASMA  CONDUCTIVITY 


Since  for  the  velocity  independent  collision  frequency 


we  may  write  (eliminating  u>) 


which  is  an  equation  of  a circle,  c - represents  the 
or  the  abscissa  (see  sketch).  Similarly, 


-(fe)  • “•  - ■>'  ■ (■£)* 


(83) 


ordinate  and 


(84) 


The  semicircle  represents  the  relation  between  the  real  and  ima- 
ginary parts  of  conductivity.  Every  point  on  the  circle  corresponds 
to  some  specific  value  of  the  operating  angular  frequency,  m.  The 
slope  of  a line  (vector)  from  the  center  to  some  point  on  a circle 
(e.g.,  Ob)  is  equal  to  — ' . It  is  obvious  from  Eq.  (30)  or  from  the 


e for  any  dependence  o: 


(86) 


This  radius  of  curvatur 
describes  not  a circle 
a displaced  circle. 


s a function  of  m,  and  thus  Eq.  (86) 
complicated  curve,  closely  resembling 


then  with  (see  Eqs.  (83)  and  (265)} 

= 4toWo)  • 

(trr)»  ♦ (5,  - ■§$)»  = (■§!)" 


(87) 
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. _ ab(ii/v)a,(ui/xi) 

? - ' 

This  normalization  makes  it  possible  for  any  functional  dependence 
of  v on  the  velocity  for  any  plasma  far  the  arc  curve  to  start  at 
3,  = 0 and  end  at  3,  = 1,  where 


For  a given  uj/v  one  can  simply  compute  5,  and  3,  from  the  above 
relation.  Results  are  shown  for  a few  cases  in  Figs.  1 and  2 . 
One  can  also  show  that  Eg.  87  can  be  written  as 

<3,  - y0)’  + (3r  - 1/2)*  = Ra 


yo  = y«(b)  = -2E>?  ijb)  . b = | , R = vVe"  + 0/4)  ■ 


If  y0  ■=  constant  this  represents  an  arc  of  a vertically  displaced 
circle  by  the  amount  y0  and  the  radius  R (see  sketch).  In  general 
this  is  a rather  complicated  curve,  with  slightly  variable  R and 


slightly  displaced  location  of  the  origin  0. 

Approximately,  we  can  choose  any  a/g,  i.e.,  y = constant  and 
draw  a semicircle  with  origin  displaced  vertically  by  y„  from  Che 


(5?  - yo)0  + Cs;  - 1/2)3  = a3 


ar  = 0,  5r  = 1 and  at  3P  corresponding  to  Che  chosen  u/v*  For 
example,  if  we  choose  b0  = wo/v  for  wh ich  5,  has  the  maximum 
(Fang38’  38a),  then 

(a,)..,  = 3, (ho)  , 

».-»<  W)  - ^ & 


0 = 9(ho)  = arc  tan  2%  . 


We  have  recomputed  oP,  o, , y0l  etc., 
2 for  integral  values  of  n and  also 


and  list  them  in  Figs.  1 and 
for  various  intermediate  n's. 


A different  type  of  plot  of  5r  • m/v  vs.  m for  n = 1 and  for 
a very  approximate  v(v)  for  air  is  given  by  Hargenau?8b 


9.  THE  ORDINARY  ELECTRIC  CONDUCTIVITY  IN  THE 
PRESENCE  OF  A STATIC  MAGNETIC  FIELD 

In  the  presence  of  an  externally  imposed  static  and  uniform 
magnetic  field  i?1  in  a plasma  and  assuming  at  first  veloaity  inde- 
pendent collision  frequency  v (weak  field  theory  is  here  used  again) 

(considering  ions  as  fixed) 

me  0 + me  J v = eE  * ftv  * H°)  ■ (88) 

magnetic  field.  The  justification  for  using  this  equation  on 
theoretical  grounds  is  only  now  becoming  apparent,  when  by  solving 
approximately  Boltzmann  equation  one  obtains  the  same  result  for 
the  case  of  velocity  independent  collision  frequency  for  weakly 

city  dependent  collision  frequency. 


f M , (89) 
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where  H2°  is  in  Gauss.  The  frequency  and  a^e  are  called  the 

electron  cyclotron  frequency  and  electron  cyclotron  angular  fre- 
quency, respectively.  Often  it  is  also  called  the  gyro-frequency. 
For  an  earth  magnetic  field  of  0.2  to  0.5  Gauss,  3cl06  < gj  • < 

IS  > 111  • for  an  earth  field  is  satisfied  for  waves  of  l < 20m,  i.e., 
f > 1.5  Mc/s.)  Then  Eq.  88  written  out  in  the  components  becomes 


(90) 


Assuming  periodic  solutions  of  the  type 


we  obtain  from  the  first  two  equations 


The  solution  of  these  equations  gives 


(91) 


(92) 

(93) 


(94) 


(95) 


(104) 


98-100  become 


(oxx'  R)  Ex  + (oxy  + B)Ey  + °xzEs  = 0 U05) 

(oyx‘  B)Ex  + (oyy  ' R}Ey  + °ysEz  = 0 (106) 

°zxEx  + °zyEy  + (ozz  ' C)Ez  = 0 • (107) 


solution  for  any  E^,  Ey , E^,  the  coeffici- 


°xy  ' *°yx  " "B 
<z  ’ °ys  = °zx  = °zy  = 0 • 


give  for  the  complex  conductivity  tensor. 


(108) 

(109) 

(110) 


(114) 


cub) 


and  thus  (6m  = 0 for  k i l and  6^  = 1 for  k = f),  using  previous 

(120) 

(121) 

(122) 

where  we  define  for  brevity  0+  = u>  + » • and  o_  = i»  - »h'  • 


y references  one  finds  only  the  real  part  of  dielectric 
) and  the  real  part  of  electric  conductivity  (a  ),  we 


(r)  _ 


(128) 


For  the  imaginary  part  of  dielectric  constant  we  obtain 

(129) 

(130) 

(131) 

(132) 

(133) 

where  we  again  write  the  connecting  relations  (k,t  = x,y,z)  as 


(r)  = -a  <r>  = . 


"'kl  - °kltr>  +*Yt(i)  = °k/r)  + "<«k/-  (139) 

_I±£  = \tr)  = «Ki  - = «k/r)*  1^(135) 

completeness,  the  imaginary  part  of  electric  conductivity 


from  Eg.  (119) 


THE  ELECTRIC  CONDUCTIVITY  OF  A WEAKLY  IONIZED  PLASMA 
THE  PRESENCE  OF  AN  EXTERNAL  STATIC  MAGNETIC  FIELD 


One  can  generalize  Che  theory18’20  of  electric  conductivity 
for  weakly  ionized  plasma  in  the  presence  of  a static  magnetic 
field  by  allowing  the  electron  collision  frequency  to  be 

compare  with  Eq.  (34)  and  Eqs.  (145)  - (149)),  (see  Appendix  H for 


( 142) 


and  where  on  Che  right  hand  sides  we  used  Maxwellian  electron 
velocity  distribution.  If  m,*  = 0,  0+  = 0 = to  and  cy  = S+  = S"  = 
S0  = o'  (given  in  Eq.  (34);.  Since  these  quantities  are  identical 
in  the  form  we  see  chat  if  we  define  o1  and  Oj  for  S0,  we  can  define 
in  an  analogous  way  a*  and  a3*  for  S+  and  o1_  and  a3~  for  S’  (see 


Eq.  (143), 


(1*8) 


From  these  expressions  using  Eq.  (119) 
a = ‘ 1ki) 

we  obtain 


°'xx  = °'yy  = (S+  + s‘> 


(149) 


(150) 


(151) 

(152) 


n show  Che  following  equalities 


£UL-(?L 


(161) 


In  Table  1-A, 


f u,/v,  where  v is  defined  by  Eq.  (61). 


Q.V  , v^wNV  . 

sfi  a a i a 

a ill  “11:11 

ad  a*  a "1 

a al  ill  s 


as  a_  a a. 

a al  II  ill 


0.8873. 


Thus  we  found  the  following  quantities: 


obtain 


0l"  = 1.6575 


= 1.2660. 
o clarify  the  u: 


A and  the  footnote). 


Example  of  Special 


°xy  " °yx  0 ' 

2.  If  tu  = 0,  <t)„  to  (D.C.  Case  with  a magnetic  field) 


fl+  = 2<u 


°*y  = °y*  = 0 ' 


approximately 
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Comparison  with  Other 
Shharofsky25 


Notations 


defined  six 


quantities;  g , h , g , h , 


Comparing  with  our  notation  one  can  show  that 


-,^o  +_>Fo 

* °5  h_  ' °3  ^hQ 

. 

Ginsburg6  uses  six  other  quantities,  K = K (n+/v),  K ( , 
K Cni/v),  Ko(0+/v),  K^( Cl_/v)  and  K^Coi/u)  defined  by  gKe(o+/v)  means 
K as  a function  of  O./u,  etc.] 


czz(l>  - 1 - P 5 K^Cui/u)  (164c  ) 


15.  ELECTRON ’NEUTRAL  EIASTIC  CRBSS  SECTIONS 


elastically  from  the  i-th  neutral  atom  or  molecule  through  an  angle 
e and  with  relative  velocity  v is  o(l)(v,e),  then  the  scattering 


cross  section  for  the  momentum  transfer,  (sometimes  referred  to  as 
the  diffusion  cross  section),  in  the  center  of  the  mass  system, is28 


directly 


(166) 
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Theoretically  one  usually  obtains  phase  shifts,  , from 
quantum  mechanics . Then 


5 Y <m  + W s 


m is  the  reduced  m. 
divided  by  2tt.  Ob 


n also  show  that 


w energies  only  a few  phase  shifts  (i.e.,  mostly  with  i = 0 
ve) , sometimes  also  with  1 = 1 (p-wave),  and  1=2  (d-wave) 
d rarely  higher)  contribute  to  the  cross  sections.*  Usually 


(i)  . 


reported  from  experimental  scattering*  data  ( e .g . , using 
molecular  beams)  as  Q = Q(^ev) , where  Egv  is  the  energy  of  electrons 
in  electron  volts.  For  electron-atom  or  electron-molecule  cross 
sections  the  conversion  relation  (for  nonrelativistic  energies)  is 

Often  the  cross-sections  are  reported  in  na  , where 


e that  if  1 = 0 only,  Q£  ' 


Frequently  experimental  data  are  given  for  collision  proba- 
bility Pc,  i.e.,  fraction  of  electrons  scattered  out  of  a beam  per  1 
centimeter  of  path  in  a gas  pressure  of  1 millimeter  of  mercury  at  0°C, 

Q(in  cm2)  = 0.02823  x 10'15  • P0  tin  ^), 


Q (in  naQ2)  = 0.322  * 
Also  1 Rydberg  = 1 Hartree  = 13.605 


= 13.605(ka  )2 


13.605ERydbergs, 


respectively.  Further, 


kao  = 0.2711  ^ 


The  approximate  range  of  energies  at  which  cross-sections 
are  needed  (for  a given  temperature)  for  calculation  of  collision 
frequencies  are  given  in  Table  5.  Thus,  we  see  that  for  T = 300°K 
we  need  mostly  cross-sections  at  very  low  energies,  i.e.,  at 
7x10  < E j 0.3  ev  and  at  6000°K  we  need  cross-sections  at  higher 

energies  0.1  < E <,6.5  ev.  Actually  one  needs  in  evaluating  the 
integral  for  collision  frequency,  also  cross  sections  at  lower  and 
higher  energies  than  indicated,  however,  the  most  significant 
contribution  to  the  integral  will  come  from  this  range  and  from 
nearby  energies  at  which  cross  sections 


have  resonances. 


Often  one  uses 
Ref.  60),  then 


13.  ELECTRON  NEUTRAL  CROSS-SECTIONS 


POLAR  MOLECULES 
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THE  E1ECTRIC  DIPOLE  MOMENTS  OF  SELECTED 
MOLECULES  IN  A GASEOUS  STATE* 


NaCl 


CsCl 


W°2 


He  (Debyes) 


0.0  - 0.16 
0.10 


C-Cl 


C=0  C-0 


, Debyes  0.4  1.6 


The  accuracy  of  the  so  determined  electron  molecule  cross- 
section  for  momentum  transfer  depends  on  the  validity  of  the  Born 
approximation. 

Since  the  derived  ~ E'1  or  ~ 1/v2  and  the  electron  neutral 
collision  frequency  v — 1/v,  we  can  use  Table  A-l  for  n = -1 
(s  = -1)  and  for  pure  gas  consisting  of  a single  species  of  polar 
molecules  we  obtain  the  effective  corrections  to  the  electric 
conductivity,  i.e.,  0^,  and  »3>  vK  = 0.8836  v and  from  Table  B-l 
for  n = -1  (s  = -1)  one  can  obtain  the  contribution  of  a given 
i-th  polar  molecule  to  the  total  electron  neutral  collision 
frequency  vg  n(w  >>  v),  i.e.,  from  Eq.  308 


(169c) 


q£|}  = 4.94  x 10'16  ue2  . (169d) 

Thus,  for  example,  for  a weakly  ionized  H„0  vapor  at  300°K  the 
electron-neutral  collision  frequency  ( IU  » 9)  is 

/u  n\  q lfi  9 (H-0) 

\ D2'  = 2-7752x10  (4.94x10  ib)  (1.82)^  N *3 

= 4.54X10'7  N<«2°>  , sec'1  . 


14.  ELECTRON  COLLISION  FREQUENCY  WITH 
NEUTRAL  ATOMS  AND  MOLECULES 

The  electron  neutral  collision  frequency,  i.e., 

the  number  of  collisions  per  second  an  electron)  makes  with  i-th 
type  neutral  atoms  or  molecules,  is* 


V,n  Cv)  = NiQm('°  ’v  • 

where  N,  is  the  number  of  i-th  neutral  atoms  or  molecules  per 
cubic  centimeter,  is  the  momentum  transfer  elastic  cross 

section  (in  cm  ) for  electron  with  i-th  atom  or  molecule  and  v 
is  the  electron  velocity  in  cm/ sec. 

The  total  collision  frequency  (for  momentum  transfer)  of 
an  electron  with  neutral  atoms  and  molecules  is 


> ■ E»a“ 


(170) 


In  principle  for  a gas  (at  high  temperature  and  up  to  a moderate  density) 
consisting  of  atoms  and  molecules  in  various  excited  states,  one 
should  consider  each  state  of  an  atom  and  molecule  as  a separate 
i-th  species,  with  the  corresponding  electron  elastic  cross 
section,  0^^  . However,  at  present  we  do  not  know  much  about 
momentum  transfer  elastic  cross  sections  for  electrons  with  atoms 
and  molecules  in  various  excited  states . 


We  are  neglecting  inelastic  collisions,  i.e.,  resulting  i 
the  ionisation , recombination,  excitation  ana  radiation. 


For  [d  » v i the  average  electron  neutral  collision  fre- 
quency is  (see  Eq.  ( 49  ))  for  electrons  with  i species  for  a 


It  is  important  to  point  out  that  not  the  collision  frequency  v , 
but  v ■ v2  must  be  averaged  over  velocity  distribution  of  elec- 
trons (R.g. , Maxwellian) . This  often  is  a source  of  a considerable 
error  and  an  origin  of  an  unnecessary  variety  of  "effective11 
collision  frequencies  (Eq.  (185), etc.).  The  case  of  v (called  by 
Ginzburg6,  Veff)  for  , » » is  to  be  preferred  since  it  forces 
the  corrections  0l,  a2  ani  o3  in  the  electric  conductivity  to  be 
unity  for  tu  » v»  and  v can  be  compared  directly  with  the  measured 

Collision  Frequency  for  Rigid  (Elastic.  Spherical)  Neutral  Atoms 
and  Molecules* 

For  an  electron  scattered  from  a neutral  particle  (N 


(171) 


particles  in  1 cnr  of  a weakly  ionized  plasma)  having  radius 
(collision  distance)  from  the  kinetic  theory  of  gases 


(172) 


ve  n = N •CiWoa)*v  . (173) 

Averaging  over  Maxwellian  distribution  (with  constant  nr  ; 

<ve>1?=  «nr02)  0,  v • (174) 

For  the  case  of  x v , from  Eq.  ( 49)  for  Maxwellian  distribution 
for  a weakly  ionized  plasma,  we  obtain  the  proper  average  (see 
next  section) 

<“> 

This  is  the  same  expression  as  one  would  obtain  from  the  kinetic 
theory  for  rigid  spheres  and  which  was  often  used  for  the  ionos- 

If  q is  the  electron-neutral  cross-section  for  the  momen- 
tum transfer  and  if  at  low  energies  Q~Q0  - constant,  then  from 
Eq.  (49) 


(176) 


Various  Definitions  of  Average  Collision  Frequency 


Often  in  the  literature  is  used  (for  a general  velocity 
dependent  cross  section}  it  has  no  significance)  the  collision 
frequency  taken  at  the  energy  equal  to  kTe,  i.e. , at 

E*  = kTe=  -i  mev*2  ; v(E*)  - v , 

at  which  the  electron  neutral  elastic  cross-section  for  the 
momentum  transfer  has  the  value 


= Q(E*)  , (177) 

For  a simple  energy  dependence  of  the  cross-section  of  the  j-th 
neutral  species  given  by 


From  Eq.  (178),  Eq.  (49) 

v =—  T(s+3)  v , 

V" 

where  the  gamma  function 


(178) 

(178a) 


(180) 


r(lM-l)  = 


( 


wVw 


(181) 


Often  we  will  use  the  following  values  of  the  gamma  function33 


r(m-l)  = m!  , for 


'GK  ' ' 

r(H*"  >'({).  if  ^ j a“> 

n«j)  ■ . 

Rnother  average  collision  frequency  often  used,  and  again 
of  no  significance  for  the  general  velocity  dependent  cross-section, 
is  the  collision  frequency  averaged  over  Maxwellian  distribution, 


instead  of  the  correct  average  v (see  Eq.  61). 

For  m » v only  P can  be  directly  compared  with  experimental 
frequency  for  weakly  and  for  fully  ionised  plasma.  (For  >jj  « v, 
only  Vj£  can  be  uniquely  defined  and  experimentally  deduced,  see 
Eq.  51).  Only  for  a simple  energy  dependence  of  the  cross-section 
v*  and  <v>  are  different  by  a constant  factor  from  v (or  , i.e., 


and  from  Eq.  (180)  and  Eq.  (184)  (omitting  j-th  superscript)  we 


(185a) 
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Q 


r3/!,a 


3^  ' 


Q 


(185b) 

(185c) 

( 185d) 
(185e) 
( 185f ) 
(185g) 


The  ratio  of  to  v for  various  n is  given  in  Table  3 


15.  CORRECTIONS  TO  ELECTRIC  CONDUCTIVITY  FOR  ft 
ONE  NEUTRAL  COMPONENT  WEAKLY  IONIZED 
PlASMft  WITH  GIVEN  Q = t^v"-1 


transfer  for  a given  one  component  gas  be  given  by 


(186) 


where  is  a constant.  This  dependence  could  be  obtained  from 

since  the  exponential  weighting  factor  in  Eq.  (35)  gives  only 
significance  to  the  low  velocity  contributions. 


ve,n  = = l*l0Vn  = NC^  (— w"/2  , (187) 

and  for  is  » v,  the  average  effective  collision  frequency  for 
a weakly  ionized  plasma  is  from  Eq.  (61)  for  Maxwellian 


(189) 
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w-n/2 


V * C'2  = T * P'1  = i 

— , cz  - s , c3=^ 


(190) 


(191) 


. (46)  and  Eq.  (47)  we  obtain  for  Maxwellian  distribution 


(192) 


/ 


W3/2e'WdW 

W"  ♦ cn!J  • 


(193) 


from  Eq.  (54)  through  Eq.  (58). 

For  integral  values  of  n (n  = -3,  -2,  -1,  1,  2,  3)  and 
and  b = l>!” — one  can  use  functions  discussed  in  detail  and 

tabulated  by  Dingle  et  al.34  For  the  case  of  n = ± 2 I1  and  Ij 
are  related  to  the  sine  and  cosine  integrals  (or  Fresnel  integrals) 
and  for  n = ± 1 to  the  generalized  exponential  integrals.  For  the 
case  of  n = ± 1 and  n = 2,  Molmud36  rederived  the  same  relations. 
Using  Dingle  et  al.  tables,  Shkarofsky  ^ tabulated  functions  g 
and  h (see  Eqs.  (72)  and  (73))  for  a few  special  values  of  b. 

suits  for  n = -3,  -2.5,  -2,  -1.5,  -1.0,  -0.5,  1,  1.5,  2,  2.5  and 
n = 3 obtained  by  numerical  integration.  For  the  integral  values 
of  n,  our  solutions  agree  with  the  values  obtained  using  Dingle 

In  the  next  section  we  will  discuss  in  detail  the  cases  of 
integral  values  of  n,  the  connection  to  Dingle  et  al.  functions  and 
the  asymptotic  relations.  Although  a direct  numerical  integration 
of  Eq.  (192)  and  Eq.  (193)  is  very  simple  and  straightforward  (see 
Eq.  (268))  and  which  we  used  and  recommend  for  any  other  computation, 
we  are  compelled  to  list  in  the  next  section  the  relations  for  the 
integral  values  of  n because  of  a large  number  of  misprints  which 
appeared  in  the  literature  while  quoting  these  relations6’14’38  and 
insufficient  number  of  terms  in  the  asymptotic  formulae  to 
guarantee  convergence  even  to  the  first  digit.38 


AJLV 


-=-HV- 


Case  of  n = -2,  Q = Qnv~ 


V-.  ff! 


<« 

%,n  * I Nt*o  ' ^ • (201) 

C£  34 

The“  (x)  function  is  related  to  the  Fresnel  integrals 

“"'VS 

’ ° ’ 0 ( 202) 

where  x = ^ u . Abre via ting  these  through  the  sine  and  cosine 
integrals 

Cii(x)  =J^[i  - C(x)] 

„ (203) 

=h«»>-V¥B-«]  ■ 


x - Ci4  (: 
x + Cij  (; 


(204) 


and  from  the  recurrence  relation 


> ”35 

and  (needed  for  n = 2 ease) 


-it1-  w]  ■ <”> 

From  Dingle*4  the  series  expansion  gives  (for  x < 1) 

1 1 
V&^(x)  = 1 1 + 4x2  - ^|  x4  + . . . | + xJZ  sin(x-t^J)  (208) 

Ct  7(X)  = ^ 1 1 - | X2  + X4  - ...  | - x2  sin(x-t^)  (209) 
and  for  x > 1, 

(I  v,  = ij  | 1 - (P^S*2)  + (^X(H-2j(pt3)(pt4)  ...  I (210) 


m Dirgle  et  al.  for  b 


£j(x)  “ ^3  (l  ' ^5  ) 


(211) 


(2X2) 
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From  Eqs.  (199),  (200), (58),  (211)  and  (212), 

Oja  | b2  (213) 

(i  - 15  1,2 ) • b ■« 1 (2M> 

»j  = ^ (l  + |i  b2  ^ , b « 1 (215) 

and  for  b » 1,  x = « 1 (from  Eq. (200) and  Eq. (209)) 


hV¥  1 


(216) 


(217) 


(219) 


e of  n = -1.  Q - QQV  ■ v = HC^v*-1 

-ttf ■%&) 


v? 

34  W 

Dingle  eC  al.  related  the  function  to  the  exponential  integral 
ST(x)  = e'E  .(x)  , (221) 


Vw  =J  ac  • 

im  Bierens  de  Haat?6  and  Batman37  for  x < 1 if  p 


% 


Fi  + TFIXF2T  - + -+1^T} 

+ eX  [-Ei(-x)]  , 

-Ei(-x)  e ei(x)  = f ^ dt  = dt  = E^x)  , 

EX(X)  = -lnvx  + X - + -^jy  - -^y 


v = 0.5772157  is  Che  Euler's  e 


«T 


' FI  + lP-l)Cp-2)  ' 
, Dingle  quotes 


(x)  = — |l  + P + PtP-2*)  + P(p2-8px ♦ 6x2) 
’ X ( P**r  (P+xT  (p+x)6 


: + ...) 


Therefore,  for  b 


#RS-i)‘2] 


»’2(x)  = -jd-x) 

,,00  - ! tl  - |»  - !#!>«  « 


(230) 

(231) 

(232) 


(233) 


The  Case  of 


-=-!<VL 


n = 0,  q = V^v. 


which  is  the  ordinary  conductivity. 

The  Case  of  n = 1,  Q = = constant,  v = NQ0v. 

This  is  a continuation  of  the  case  discussed 
(172)  through  (176). 


/ "« ?»v ' * 


if 


e function  was  described  previously  in  the  n = 

2f 


(234) 

(235) 

(236) 


(237) 

(238) 

(239) 

(240) 


(241) 


CORRECTIONS  TO  THE  ELECTRIC  CONDUCTIVITY  FOR 
RIGID  SPHERICAL  ATOMS  AND  MOLECULES  IN  A 
WEAKLY  IONIZED  PLASMA* 


for  n = 1. 


# 

a 


This  case  was  investigated  bv  Fang 
relation  to  Dingle  integrals^ 


« * 


(246) 
who  first  pointed 
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The  Case  of  n = 2.  Q = Q_v,  v = NQ0v2. 


«e,n  aSV%-mf  • 


The  VS*  function  was  discussed  in  n = -2  case. 


S 00  * | [l-  4 *2  +^>/x  COS  (X  + =)] 

5,,,“  * i [*  * - Vs  «J>  ft  - Ji] 


(247) 

(248) 

(249) 

(250) 

(251) 


(252) 


(Mote  an  error  in  Molmud55  in  Eq.  A-7,  7 should  read  63). 


(254) 

(255) 


" 717 


(257) 


(250) 


The  Case  of  n = 3,  Q = Q^v  , v = NQQvJ 

• <*» 

There  are  no  convenient  asymptotic  series34  for  Cf  (x)  function, 
and  series  quoted  by  Dingle34are  very  slowly  convergent  (see 
discussion  after  Eq.  (198)). 


FIGURE  3 Correction  ai  for  Weakly  Ionized  Plasma 


FIGURE  4 Correction  ota  for  Weakly  Ionized  Plasma 


FIGURE  5 Correctic 


for  Weakly  Ionized  Plasma 
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The  DC  and  Related  Limiting  Cases  of  nt  =0 

For  b = = 0,  (DC  case),  (or  b = “ * "fr*  = 0, 

see  footnote  to  Table  A-l)jrfe  have  in  the  limit  of  b - 0 


From  Eq.  (58)  we  obtain  in  the  limit  of  b - 


The  direct  current  (d.o.)  electric  conductivity  is  given  by 
(see  Eq.  (60)) 


where  v is  the  average  collision  frequency  at  , » u and  v^q  is 
the  average  effective  collision  frequency  at  tu  = 0.  Thus  (see 
Eq.  (76)) 


v is  r(~^n)  r{szny  ~ ^ = g,,°°  ^ <265) 

where  g^  is  the  notation  used  by  Shkarofsky25  and  yg  by  Spitzer^6 
(see  Eqs.  69  and  76)  for  fully  ionized  plasma. 


for  determination  of  c’s  for  |®  ± m,*|/v  near  aero.  As  we  see 


• ■»> 


v - L,  ve>n  + *e,i 
Eq.  (34),  etc.,  must  be  abandoned. 

« ■ '-k - v**”W-  £ ■■(¥) ■#)  • 

i.e.,  for  n = 3,  o°o”  = ” , = 5n  slmil„ily 
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FIGURE  6 The  Ratio  of  Effective  (DC)  Collision  Frequency  ot  to-  0 to  on  Average 
Collision  Frequency  for  ie»  v as  a Function  of  n,  where  u = Avn 
for  Weakly  Ionized  Plosma  e'n 


16.  NUMERICAL  EVALUATION  OF  INTEGRALS 
FOR  WEAKLY  IONIZED  PLASMAS 


of  evaluation  of  integrals  for  weakly  ionised  plasma  (Eqs.  (46) 


for  collision  cross-section  approximated  by  finite  series 


functions  and  exponential  integrals,  both  of  complex  arguments. 

However,  the  resulting  integrals,  even  for  multicomponen 
weakly  ionized  plasma,  can  be  very  simply  and  accurately  evalua 
numerically  by  the  Simpson  formula,  or  using  n point  Laguerre39 
formula  based  on  the  roots  of  the  Laguerre  polynomials,  i.e. , 


J f(W)e”WdW  = ^ \ f(Wk)  , (268) 

where  Ak's  and  W^'s  are  precomputed  by  Salzer  and  Zucker39  and 
listed  in  Table  4 for  a n=5  point  formula.  Table  4 illustrates 


point  formula 


92 


Te(°K) 

Ek(ev)  = jmr  “k  ■ « V " E W V 


<v  . 


. is  the  electron  ion  collision  frequency  (see  Ch.  II,  ReJ 
e Qj  is  the  momentum  transfer  elastic  cross-section  01  an 
with  jth  neutral  species  and  = N../£n.  . Then,  for  a gi 
where  0 = |is  ± <ii„'  I , from  Eqs.  (46),  (47)  and  Eq.  (49), 

(Vv)  ^ 


----H  \ 


<V»>  + 1 
:E  \ '’k  wk 


and  from  Eq. 


93 


for 


for  (Jjl  =|oi 


EXAMPLES 


GAUSS-IAGUERRE  QUADRATURE  FORMULA 


= J f(W)eWdW=^  Akf(Wk) 

A.  A. 


1.52934 

3.J233S 

1.63173 


V- 


0.50000 


15  point  to 


f(0)  = 
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On  Uncertainties  in  the  Evaluation  of  Integrals  Involving  Electron- 


The  uncertainty  in  the  electron-neutral  elastic  cross  sections 
for  the  momentum  transfer*  is  about  ± 0.05  x 10"15om2  from  a given 
experiment  (e.g.,  molecular  beam).  The  initial  collection  of  re- 
sults from  the  various  methods  of  measurements  of  cross  sections 
by  various  investigators45  and  from  the  various  methods  of  calcula- 
tions give  different  values  (in  the  range  of  0 to  2 ev),  e.g.,  a 
factor  up  to  7 for  N2,  of  S for  02,  of  2 for  NO,  of  3 for  N,  of 
9 for  0,  of  9 for  ft,  etc.  However,  after  careful  review  one  decides 

appendix  G).  Then  the  remaining  uncertainty  at  best  is  still 
± 0.05  x 10'15cm2.  (assuming  constant  with  energy  £Q  = ± 0.05  x 
10_15cm2  gives*  Av  =■  * 0.2  x 1011o/c0  at  T = 300°K.)  This  limits 
the  accuracy  of  integration  at  most  to  2 or  3 significant  digits. 

The  five  point  Laguerre  formula  (Salzer  and  Zucker39)  is  definitely 
sufficient  (see  Table  <1).  Table  5 indicates  the  values  of  energy 
E , at  which  the  cross  sections  are  needed  for  the  evaluation  of 
the  integrals  (assuming  5 point  formula  is  sufficient)  for  various 


See  Table  6,  s=0,  Y = 3.BS5X1026,  Av  = 3.85Sxl026  • Sxl0‘17p/po 


3-SECTION 


17.  CONTRIBUTION  OF  A RESONANCE  IN  A CROSS 
TO  COLLISION  FREQUENCY 

Let  Che  momentum  transfer  elastic  cross-section  for 
electron  with  neutral  atom  or  molecule  be  split  into  two  parts: 
one,  without  a resonance,  which  we  approximate  by  a polynomial 
fit  in  energy  and  a resonant  part  at  E = E , which  we  approximate 

Q(E)  = £ qsEe*  + Qre-°<E-Eo)  . (278) 


Then , for  w » v from  Eq.  61  and  from  Eq.  180 , 

’••'■'w'V?  rt-*s> ia?s> 

where  the  contribution  of  resonance  cross-section  to  the  colli- 
sion frequency  is 


FIGURE  7 Assumed  Goussian  Profile  for  fhe  Resonance  in  o Cross-Section 


F(W)  = exp  [-ot2(W  - V 
After  an  exact  integration*  we  have 
1 3 


I = e'^o  - 


where  Erf  (x)  is  the  error  function 
Erf(x) 


ot2+2(|  - otE  )2  /|-otE  \ 
f e"2  du  = 1 - E2(x)  , (284) 


)=^{S'U  ' 


If  we  define  t: 


AE  by  (see  sketch) 


«E0  - ^AE)  = Q(E0  + -|aE)  = Qre-1 


AE  = — . Thus  for  AE  s 1 ev  , i)  4 , ( 287 ) 


f * ■ 1?  *f~  “ ($)  ■ ™»> 

(This  integral  can  be  obtained  by  a substitution  u = x ^ + 3 /s/p  , 
forming  e_u2  and  using  J ue"u  du  = -i  e”®2  and  J u2e~u2  du  = 


frequency  at  the  energy  equal  to  -|  kl,  i.e.,  at  the  temperature 


T.('V)  = 4,642.«  E0(ev)  . 

E0  a 0.02  - 0.06  ev  , verjjs  has  maximum  at  T ~ 100°K  - 300°K 
Eoa  0.2  - 0.65  , ve”S  has  maximum  at  T ~ 1000  - 3000°K 


Eo-£-2  ev  * has  maJcimum  f°r  10,000°K. 


r 


(AE  < 1 ev)  (290) 
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FIGURE  9 


constant  V i.e.,  if 


! .raSL'va  ' 


3 -ertEo  = xfc/5- 
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Numerical  Fits  to  Elastic  Neutral  Cross  Sections  and  the 
mMfeution  ot  cs  ss  Section  wear  zei  Energy  to  tne  Avera  e 


The  range  of  possible  power  exponents  of  v = Avn,  i.e., 

= Q0''n"  can  be  approximately  deduced  from  the  following  con- 
sideration. Since  (Eg.  (168)) 

= ^5  E (U1)  s jjs  s » * =f  (**) 

it  follows  that  with  a finite  A = ( 1+1) 

Qv2  s ^5  A , (295) 


°cm2  • Eev  s S'2  x • (296) 


At  very  low  energies  (E  < 1 ev)  only  few  partial 
to  s,  mostly  the  s wave  (1  = o)  and  the  p wave  (i 
only  approximations  (from  zero  energy  to  same  E 
ments  are  available)  with  n = -l,  o,  1,  2,  etc. 
all  energy  range,  but  not  with  n < -1,  i.e.  not  ■ 


waves  contribute 


are  permissible  for 


V = «o 

en  though 


an  be  greater  than  any  constant 
fits  up  to  zero  velocity  are  not 


possible  generally,  one  is  justified  in  using  any  power  fit  to 
experimentally  observai  cross  sections  because  very  low  velocities 
(energies)  do  not  contribute  significantly  to  the  value  of  the 
needed  integrals  (e.g.,  in  Eq.  (35)  and  (61)  we  multiply  Q by 
vsexpC-|mv2/kT)  and  integrate  over  v. 

At  higher  energies,  contributions  from  higher  partial  waves 

fit  to  the  actual  cross-sections  up  to  high  energies,  beyond  which 
the  exponential  weighting  factor  in  Eqs.  (35)  and  (61)  would  re- 
duce any  difference  in  the  value  of  the  definite  integrals  between 


energy  leaving  the  uncertainty  in  the  range  0 s E < E,.  If  we 


% = ^ . ("  = 1-2P,  v = Av") 


0 « E <:  E1  , ('.97) 


«e,n  = N (iTeS)P  % / 1"2-Pe’W*  > <2B8) 


(300) 


(301) 


= /V* 


n = -1), 

«.<*>  = 11  «Bi>  • W1[l-(1+W,)e-W']  . (3o2) 

number  of  partial  wave  contributions  cannot  go  to  sero  faster  than 
1/E,  and  the  above  expression,  it  can  be  shown,  will  give  for 
0 < n s 1 maximum  contribution  to  the  collision  frequency  from  the 
0 to  Ej  energy  interval.  (It  is  interesting  to  note  that  only  for 
p s 5/2,  i.e.,  n z -4,  one  obtains  finite  iv  n>) 


18.  AVERAGE  ELECTRON  COLLISION  FREQUENCY  IN 
TERRESTIAL  ATMOSPHERE 


quency  for  an  upper  atmosphere  as  a function  of  an  altitude  is 


the  uncertainty  of  a composition  and  electron  temperature  at  high 


•5 


m.v3 

ea 


where  \t  is  the  Debye  length, 


3279  Wf,  ' 

•at.  * 


(«»y) 


109 


no 


in 

where  T and  T.  are  electron  and  ion  temperatures,  respectively, 
is  the  electron  concentration  per  cm3  and  where 


C CONSTANTS  OF 
WEAKLY  IONIZED  GASES  FOR  VARIOUS 
VELOCITY  DEPENDENCES  OF  ELECTRON 
NEUTRAL  COLLISION  FREQUENCY 
(MAXWELLIAN  VELOCITY  DISTRIBUTION) 
-3*  (v  = Aun,  q = = q^"1,  n = 


al 

|:S8 

a2 

ZEZ 

1:1 

5:1)1; 

S.2.S1 

4.3919 

wbs 

if 

w:;; 

SiSa 

liiiii 

liiSH 

|:j«S 

1 : 0676 

1.0146 

*9997 


If  | = b,  then  «a  = a2>  «,  = a3,  («3  = ^). 

If  = b,  ^en  «+  = «!■  «+  = a3,  («2+  = «2  = «2>- 

Ifl"  1=  b>  then  “l'  = al»  “s'  = a3>  la2  = “2  = a2> 


fl2 


■Mil 


TABLE  A-l  (Cont) 


:/v  - cri/a33  - •g  = 0.6000 


ft-1  (Cont) 


voc/^  = Oj/ffa 
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9656 

9712 

9735 

9757 
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i 

.soo 

.600 

:!!! 

| 

I;-::: 
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TABUS  A-l  CCont) 


voe/v 


TABLE 


9213 

9J1B 

922« 

9261 

93*2 

9*23 


v»c/v 


0.8836 


I 
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TABLE  A-l  (Cont) 


955? 

9720 

9810 


.7543 
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TABLE  A-l  (Cont) 


1.1912 

1.2242 

1.010/ 


6741 

7113 


9220 

9298 

9365 

9931 

9961 


123* 
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200 


'loe/«  = 


APPENDIX 


AVERAGE  ELECTRON-NEUTRAL  COLLISION  FREQUENCY 
FOR  WEAKLY  IONIZED  GASES 
(MAXWELLIAN  VELOCITY  DISTRIBUTION) 


For  a given  i-th  neutral  species  if  its  electron-neutral 
(momentum  transfer)  elastic  cross  section  dependence  on  energy  (i 
velocity)  can  be  approximated  by  0^^ 


e constants , then  one  can  use  Table  6 , 
e f or  various  n ( and  5 ) and  various  tempera- 
tures, the  following  quantities  are  listed:  Ys^*',  DS^K  Sn  and 
which  are  defined  by  the  average  (in  the  sence  of  Eq.  (61)) 
collision  frequency  of  an  electron  with  i-th  neutral  species  (to 


r for  ai 


e S) 


ec-)  = D/^V 


= s C c (i)  = R (i)c  (i>P  (i>  ■ v =Fv  (i) 

&nLit  pQ  Ks  mmHg  ’ 'e,n  e,n 


The  quantities  were  computed  from  Eq.  180. 
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- fe)“,z  t ■ ■ 


Eq.  (278)). 


■ Table  B-l,  Ds  = Dg(t)  we  have 


= N * WT> 


v. ' V'*"’ 


B.10 
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where  LQ  = 2.6872  x 1019  cm"a,  and  where  1 ev  = 11605°K,  i.e.,  if 
T is  in  °K  and  E in  electron  volts 
. m v8 

■5  = W = Eev  ' B'U 

Illustrations  for  the  Use  of  Table  B-l. Example  1, 

The  momentum  transfer  elastic  cross-section  for  electron 
with  helium  was  neasured4*  to  be* 

Q = 5.3  x 10_1°cm3  for  0 < E < 0.06  ev.  B-12 

Assuming  this  to  be  valid  for  higher  energies  we  obtain  an  average 
electron-neutral  collision  frequency  (uj  >>  v)  for  weakly  ionised 
pure  helium  at  T = 300°K  and  any  pressure  (see  Table  9). 

We  have  for  s = 0 at  T = 300°K  from  Table  6,  RQ  = «. 62x10" . 


5e  n = 4.62  « 10"  • 5.3  x lo”: 
Experimental  values  for  helium  for  P 


Example  2. 

L.C.  Chen42’45  found  that  the  experimental  electron  collision 
frequency  with  neutral  neon  atoms  at  low  energies  can  be  fitted  with 
Q1  = 2.54  « 10"'“  , cm8  B-14 


TTor  He  (Table  E-l). 


Appendix 
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Find  » v)  for  Tg  = 100,  200,  300,  400,  500  and  600°K  and 

compare  with  Chen's  experimental  values. 

From  our  TableB-lfor  n = 2 and  n = 1 (i  = i,  £ = 0)  we  fi 


T°K 

100 

200 

300 

400 

500 

E00 


Qjtenr)  = 2 
5e,n  = qV^/aP 


VTv  ■ 


T °K 


(Note  that  Chen  used  Po  from  P/P  = T/273.15.) 

Using  the  other  expression  for  the  cross-section 

Q^cm2)  = 1.07  x 10"17  + 2.17  x 10_1^E^  , 

we  find 

^e,n  + t,1/aRl/a^  PmmHg’  ‘’o  ~ i-O^lO  , <h/a=  2 


>.17x10 


FIGURE  B-l  Electron-Neutral  Collision  Frequency  for  Neon 
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the  electron-neutral  collision  frequency  for  very  high  mirco- 
wave  frequencies,  for  equilibrium  temperature  of  300°K  and 
concentration  of  neutral  molecules  of  10 19  cm"3  (electrons  are 
assumed  to  be  present  from  an  external  source).  Answer.  Using 
Table  B-l  for  S = -3/4, 


(0.03xlO"1S)  -(2X1019)  = 7. 
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APPENDIX  C 


CORRECTIONS  a1,  «z  and  a FOR  THE  ELECTRIC  CONDUCTIVITY 
FOR  WEAKLY  IONIZED  N,,  0,  AND  AIR 


1S1 
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FIGURE  C-l 
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FIGURE 
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APPENDIX 


AVERAGE  ELECTRON  COLLISION  FREQUENCIES  FOR  WEAKLY 
IONIZED  N2,  Oj,  NO,  N,  0 AND  DRY  AIR 


N(cm“* ) - l„  p/Po 


1»  = 2.6871xl019  , 


p..««  = 760(p/po)T/273.16  . 

OoCNp)  = 1.250x10“*  gm/cm*  ; pp(dry  air)  = 1. 223x00"*  gm/cm* 

Po<l%)  = 1.428x10-*  gm/cm*  ; pp(h.temp.air)  = 1.293x10'*  gm/cm* 

Dry  air:  78.0881*  N,  ; 20.9795*  0,  ; 0.9324*  A. 

High  temp.  air:  78.084*  Ns  ; 20.946*  a, ; 0.937*  A;  0.033*  CO, 
(only  concentrations  of  N, , Q, , N,  0,  NO  and  electrons  were  used). 


tabulated  in  the  next  tables  and 


5,W  = ypW/2. 
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FIGURE  D-l  Average  Collision  Frequency  for  Weokly  Ionized  Dry  Low  Temperature  Aii 
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FIGURE  D-3  Average  Electron  Collison  Frequency  for  Atomic  Oxygen  and  Atomic  Nitrogen  (al  >l 
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FIGURE  D-5  Total  Average  Electron  Collison  Frequency  of  Equilibrium  Air  (a)>>y) 
(Electron-Neutral  t-  Electron-Ion) 


AVERAGE  ELECTRON  COLLISION  FREQUENCY 
HIGH  TEMPERATURE  EQUILIBRIUM  DRY  A: 


2000 


8.459x10° 

8 . 4B4  xlO7 
8.492x10° 
8.495x10° 
B.495xl0!o 


459x10° 

492X10° 

495xlQ10 


8.459 


8.495 

8.495 


1.127x10’ 

1.133X10° 

1.148X10° 

1.173x10” 
1.176  xlO1  ° 
1.177X1013 

1.299x10° 
1.307x10° 
1.322X1010 
1.340x10” 
1.349x10' 2 
1.352X1013 

1.490  XlO10 

1^523'XlO13 
1.  53HX1013 


9. 925x10’ 
l!oO8Xl010 
1.009x10” 

1.010X10’- 3 

l.oio xio13 

1.132x10’ 

l!  150x10° 
1.165  xlO1 0 

l!  176  xlO1  2 
1.177  XlO1 3 


1.317x10° 
1.316x10° 
1.326  XlO13 


1*526  xlO1  ° 
1.531  XlO13 


front  calculations 


by  iiochstim. 


(Coned) 


S96xl0j  = 


126x10’ 

203x10° 

242X10° 
370X1010 
835x10” 
350 XlO1  ° 


113X10° 

124 XlO10 

09SX1011 


277x10” 

217X1013 


333XlOl° 
780 XlO10 
871x10” 
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TftBLE  D-l  (Contd) 


030xl07 

285x10= 

381x10= 

650x10” 
207x10*  s 
6 39x10*  3 


097  Xl0'° 

015x10** 
040X10'* 
911X10*  = 


APPENDIX 


AVERAGE  ELECTRON  COLLISION  FREQUENCIES  FOR 
WEAXLY  IONIZED  He,  Ar,  Kr  and  Xe 


Tables  (Appendix  E and  F)  only  two  significant  digits 
cted  to  be  accurate;  the  remaining  two  digits  are 
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TABLE  E-l 


AVERAGE  ELECTRON  COLLISION  FREQUENCY 
FOR  WEAKLY  IONIZED  HELIUM 


IIIIIIIIIIIIIIIIIIIIIIII 


(»»  v) 
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AVERAGE  ELECTRON  COLLISION  FREQUENCY 
FOR  WEAKLY  IONIZED  KRYPTON 


10000 

12000 
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AVERAGE  ELECTRON  COLLISION  FREQUENCY 
FOR  WEAKLY  IONIZED  XENON 


_5l_ 


13000 

15000 


18000 

20000 


17.20 
11  "99 


7!  369 

4!  804 
3.924 
3.709 

4!  550 


34.33 


1.151 

0.7052 

0.5332 

0^4672 


0.6761 

0.7826 

ol9742 

1.054 


1.226 


1.316 

1*344 

1.350 


I 


175 


5)  ,01 


APPENDIX  F 

DC  EI£CTRIC  CONDUCTIVITY 
FOR  Oj,  N2  AND  DRY  AIR  (T  s 2000°K  ) 
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V»c  = vSc  • 10”  £ > V.C  = • *<°>  = Ob  at  „ = 0,  etc. 
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APPENDIX 


REVIEW  OF  ElASTIC  CROSS  SECTIONS 


history,  both  in  the  value  and  in  the  shape,  which  only  after  pro- 
longed and  repeated  measurements  by  different  methods  and  detailed 


References  are  at  the 


this  Appendix. 
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Refs.  1,  2,  and  Bruche,  Ref.  8) 
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190 
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fical  integration)  of 


Comments  on  the  Adapted  Values  (for  the  nuraer 
Electron  Momentum  Transfer  Cross  Sections. 


From  review  by  Shkarofsky,  Bachynski,  and 
Johnson1  (1961). 

[Based  on  private  communication  to  Shkarofsky 
et_al.  by  R.W.  Crompton  and  L.G.H.  Huxley 
(subject  to  a correct  distribution)  using 
data  of  R.A.  Nielsen  and  N.E.  Bradbury,  Phys. 
Rev.  51,  69  (1937).]  The  numbers  were  read 
from  Fig.  2a  in  Ref.  1 and  extrapolated  to 
energies  below  0.025  ev. 

Bruche8  (1927).  Shape  verified  by  Neynaber 
_et_al. 14  (1961). 


3lpsJ  (1961) 
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io-“  .nv7 


this  gives  Q = 0.426x10'“  cm3  at  E = 0.05  ev 
Q = 0.129x10-“  at  E = 0.003  ev. 
By  taking  straight  line  approximation,  so 
that  its  slope  is  equal  to  the  slope  given 


(1961)],  who  used  data  of  R.A.  Nielsen,  Phys. 
Rev.  50,  950  (1936).  Read  from  graph  in 


1964) 


slight  departure  in  our  Fig.  1,  which  has 
no  noticeable  effect  on  the  collision  fre- 
quency). 


(Actually  they  measured  the  ratio  of 
Q(0)/Q(a> ) and  used  values  for  Q(q,)  of 
E.  Briiche9  (1927), 


Neynaber,  Marino,  Rothe  and  Trujilo' 
(1963)  [actually  measured  ratio  of 
Q(N)/Q(N„)  and  used  Q(Na)  from  Normand11- 
Bauer  and  Browne  8 (1963) 

Joined  smoothly  between  Refs.  7 and  8. 


Based  on  Bailey  and  Sommerville  15  (1934). 


Extrapolated  from  above. 

Read  from  graph  in  Massey  and  Burhop  5 
(1956).  Based  on  data  obtained  by  Ramsa 
method. 

Joined  smoothly  between  values  from  Refs. 


and  Phelps  (1961),  and  Engelhardt 
Phelps  ” (1964).  Argon  was  used  for 
for  T s 2000°  only. 


Pack 


TABLE  G-l 

ADAPTED  FOR  CALCUIATIONS 

ELECTRON-NEUTRAL  COLLISION  CROSS  SECTIONS  FOR  MOMENTUM 
TRANSFER  AS  A FUNCTION  OF  ELECTRON  ENERGY 

(Q  in  units  of  10-1Qcma) 


E(ev)  Q(0a)  Q(Ng ) «NO)  0(N)  Q(0) 


.2280 

.2410 

.2530 

0.2640 

.2780 

.3080 

0.3212 

.3320 


1.204 

1.215 

1*230 
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.3690 

.3770 

.3840 

.*3980 

.4020 
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E(ev) 


Bgssfegssssslsssgjiasalsisi ssesssslsss 


QC  NO  ) 


609 

610 
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5- 
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rsi,r“" 
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15.  Bailey,  V.S.  and  J.M.  Somerville,  Phil.  Mag.  17,  11G9  (1934; 
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